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In this paper, we use momentum sum rule to compute the frac- 
tions of momentum carried by quarks and gluons in a self-similarity 
based model of proton. Comparison of the results with the prediction 

of QCD asymptotics is also made. 
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1 Introduction 

Self-similarity is a possible feature of multi-partons inside a proton at small x 

^ suggested by Lastovicka pQ and pursued by us later (2J [3] . We use momentum 

i—< sum rule [I] to compute the fractions of momentum carried by quarks and 

>• gluons in such a model. As this sum rule needs information about x in the 

regime < x < 1, we extrapolate the validity of the model beyond its original 

range 6.2 X10 -7 < x < 0.01. In Section 2, we outline the formalism and in 

Section 3, we discuss the results. 



2 Formalism 

2.1 Self-similarity based model at small x 

The model is described by self-similarity based parton density functions 
(PDFs) in Ref. [1]: 
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where the parameters are 



D 1 


= 0.073 ±0.001 


D 2 


= 1.013 ±0.01 


D 3 


= -1.287 ±0.01 


Ql 


= 0.062 ± 0.01 GeV 2 
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i.e. 



qi (x,Q 2 )=e D »f(x,Q 2 ) 
where D l Q is the only flavor dependent parameter and 
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is the flavor independent function for describing x and Q 2 dependence of the 
PDF. 
Defining structure function 



F 2 (x, Q 2 )=xY, e 2 {q t (x, Q 2 ) + q t (x, Q 2 )) 
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we get 
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i.e. 



F 2 (x,Q 2 )=e D °xf(x,Q 2 ) 
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(6) 

(7) 



where 

n f 

with D = 0.339 ±0.145. 

2.2 Momentum sum rule 

The momentum sum rule [I] is given as 
l 



f x^( qi (x,Q 2 )+qi(x,Q 2 )) dx + I G(x,Q 2 ) dx = 1 (9) 





where 

i 



(*>«= /X)(9i(^Q 2 )H-9i(^Q 2 )) <& (10) 



o 
and 
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are the momentum fractions carried by quarks and gluons respectively. 
Using Eq (3) in Eq (10), 



(x) q = e 3 ° I xf(x,Q 2 ) dx (12) 

o 

where 

n f 

e^ = ]jr(e D ° + e D °) (13) 

i=i 

Note that D and D are not identical except in specific model [51 El E] , where 
partons have integral charges. Hence the momentum sum rule as given in Eq 

(9) becomes 

i 

[ (aF 2 (x,Q 2 )+G(x,Q 2 )) dx = I (14) 



where 
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a =- ^ ( 15 ) 
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which has to be determined from data. 

Since D 3 is negative and D x is positive (Eq (2)), there exists a singular- 
ity in x at x ~ 0.019 in the structure function F 2 (x,Q 2 ) (Eq (6)). Any 
viable structure function is not expected to have such a singularity in the 
region < x < 1. We therefore evaluate the integral of Eq (14) by Cauchy's 
principle value integration [8] leading to 

(x) q = a.I (x min ,Q 2 ) (16) 

where I (x min , Q 2 ) is an infinite series. 
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where 



a = D 2 -l + D 1 \og(l + % 

(3 = D 2 -l 

H = i- a = l- (D 2 -l + D 1 \og(l + ^ 

p = l-p = l-(D 2 -l) 
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Here x min is introduced to take care of the end-point singularity x = in the 
model. Taking only the first term of the integral I (x m i n , Q 2 ), we have 

(x) q = o.(0.024607).log[ ymaa; ~ yo [.^ (19) 



3 Results 

In order to make numerical analysis, we have to determine a and x m i n defined 
in Eq (15) and Eq (19) respectively. We take 

l 

(x) q = 0.465 ± 0.023 and / F 2 (x, Q 2 ) dx = 0.148 ± 0.007 (20) 

o 

at Q 2 = 15 GeV 2 [H] and evaluate 



( x )i 



J F 2 (x, Q 2 ) dx 



3.1418 



(21) 



Using it in Eq (19), we obtain 
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3.4953 x 10 



-4 



(22) 



In Table 1, we record the values of (x) q and (x) g for several representative 
values of Q 2 . We also show the values with a = 1 corresponding to partons 
having integral charges ±1 in order to explicitly show the effects of charges 
in the partition of momentum among quarks and gluons in the model. 

Table 1: Values of (x) q and (x) g for a few representative values of Q 2 . 



Q 2 (GeV 2 ) 


a = 1 


a = 3.1418 


a = 1 


( x )g 
a = 3.1418 


10 


0.0716 


0.22516 


0.9283 


0.77484 


20 


0.1433 


0.45032 


0.8567 


0.54968 


30 


0.2149 


0.67548 


0.7850 


0.32452 


40 


0.2867 


0.90064 


0.7133 


0.09936 


45 


0.3225 


1 


0.6775 





50 


0.3583 




0.6418 




60 


0.4299 




0.5700 




65 


0.4658 




0.5342 




70 


0.5017 




0.4983 




75 


0.5375 


> 1 


0.4625 


<0 


80 


0.5733 




0.4267 




90 


0.6449 




0.3550 
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Figure 1: (x) g and (x) g vsQ 2 (for a = 1). 
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Figure 2: (a;) (? and (x) g vsQ 2 (for a = 3.1418). 

In Figure 1 and Figure 2, it is shown that within the leading term ap- 
proximation used, (x) q increases with Q 2 while (x) g decreases. For a = 1 
(Figure 1), (x) q « (x) g at Q 2 w 70 GeV 2 and for a = 3.1418 (Figure 2), 
(x) q ~ (x) g at Q 2 ~ 22 GeV 2 i.e. both quarks and gluons share momentum 
equally. For a = 3.1418, the pattern is the same but the share of momentum 

for quarks increases threefold. It is to be noted that one usually expects 

i 

the other pattern, i.e. (x) g = jG(x,Q 2 ) dx should increase with Q 2 and 



( x )<? = J-^2 (x,Q 2 ) dx should decrease [10]. This feature is presumably due 

o 
to the overestimation of the large x quarks in Eq (6) and the crude ap- 
proximation of taking only the first term of the infinite series as in Eq (19). 
Further, as Q 2 increases, (x) q ps 1 at the saturation scale of Q 2 ~ 140 GeV 2 
(for a = 1) and Q 2 ~ 45 GeV 2 (for a = 3.1418) respectively, beyond which 
the momentum sum rule breaks down. Thus the model in the original ver- 
sion has inherent limitation because of the breakdown of the momentum sum 
rule beyond the regime 45 < Q 2 < 140 GeV 2 . For comparison, we record the 
experimental values of (x) q and (x) g in Table 2 available in literature. 



Table 2: 



(X)q/(x)g 


Value 


Q 2 (GeV 2 ) 


Reference 


(X)g 


0.3 ±0.008 


3.5 


[10] 


(X)q 


0.465 ±0.023 


15 


[9] 


\X)q = \X)g 


0.50 


10-20 


[11] 


(*)g 


0.43 ±0.03 


30 - 100 


[12] 



We also compare the pattern of change of (x) q and (x) g of the present model 
with standard QCD. Asymptotically (Q 2 — > oo), perturbative QCD predic- 
tions of (x) q and (x) g [13j [H] are 
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(24) 



where nj is the number of active flavors and n g is the number of gluons. For 
SU(3) C , n g = 8. 

In Table 3, we evaluate (x) q and (x) g for various flavors rif = 3, 4, 5, 6 setting 
n g = 8 using Eq (23) and Eq (24). 

Table 3: Flavor dependence of asymptotic prediction of (x) q and (x) g 
in perturbative QCD. 



n f 


( x )g 


( x )g 


3 


9/25 


16/25 


4 


3/7 


4/7 


5 


15/31 


16/31 


6 


9/17 


8/17 



This table shows that only for rif = 6, (x) q > (x) t 



In the present model, on the other hand, such feature sets in even at 
much lower Q 2 range as discussed earlier. The study of flavor dependence 
of (x) q and (x) g in the present model is not possible as it is absorbed in the 
model parameter D Q (Eq (8)). The effects of large x and higher order terms 
are currently under study. 



It is also interesting to note that the asymptotic partitions of quark and 
gluon momenta of Eq (23) and Eq (24) have been questioned in recent liter- 
ature [T3]. The suggested alternative asymptotic partition is 

lim (x) q = f (25) 

q 2 ^oo n g + Grif 

lim (x) g = - * (26) 
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leading to (x) g ~ —{x) q and (x) g ~ 7\( x )<i ^ or n f = ^ anc ^ n f = ® respec- 

-Lo y 

tively even asymptotically. In the present model such specific partitions set 

in at finite Q 2 . Thus the main feature of the self-similarity based present 

analysis, i.e. (x) q > (x) g with increased Q 2 might conform to the QCD 

expectation as well if the result of Ref. [15] is taken seriously. 

The model, even though predicts such interesting results, has several de- 
ficiencies of its own. First, the model has phenomenological validity only in a 
limited x-range. Second, the model falls short of accommodating momentum 
sum rule for large Q 2 , Q 2 > 140 GeV 2 . Finally, as it is purely self-similarity 
inspired, the model does not conform to standard DGLAP evolution. 



4 Conclusions 

In this paper, we have used momentum sum rule to compute the fractions 
of momentum carried by quarks and gluons in a self-similarity based model 
of proton. The main feature of the present analysis is that the momentum 
fraction of quarks increases with Q 2 , while the corresponding fraction of glu- 
ons decreases. This prediction conforms to asymptotic QCD expectations, 
as reported recently [TS]. The difference is that such momentum partitions 
set in at finite Q 2 in the present model. 
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